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Abstract 

y—i ■ 

The polarization of the Cosmic Microwave Background (CMB)is a powerful obser- 
ly-j | vational tool at hand for modern cosmology. It allows to break the degeneracy of 

fundamental cosmological parameters one cannot obtain using only anisotropy data 
and provides new insight into conditions existing in the very early Universe. Many 
£0 [ experiments are now in progress whose aim is detecting anisotropy and polariza- 

tion of the CMB. Measurements of the CMB polarization are however hampered 
by the presence of polarized foregrounds, above all the synchrotron emission of our 
S^-f Galaxy, whose importance increases as frequency decreases and dominates the po- 

larized diffuse radiation at frequencies below ~ 50 GHz. In the past the separation 
-£3 . of CMB and synchrotron was made combining observations of the same area of sky 

at different frequencies. In this paper we show that the statistical properties of the 
polarized components of the synchrotron and dust foregrounds are different from 
the statistical properties of the polarized component of the CMB, therefore one can 
build a statistical estimator which allows to extract the polarized component of 
the CMB from single frequency data also when the polarized CMB signal is just a 
fraction of the total polarized signal. Our estimator improves the signal/noise ratio 
for the polarized component of the CMB and reduces from ~50 GHz to ~20 GHz 
the frequency above which the polarized component of the CMB can be extracted 
from single frequency maps of the diffuse radiation. 
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1 Introduction 



Almost a decade elapsed since the first detection of the anisotropy of the Cos- 
mic Microwave Background at large angular scales (> 10°) [Strukov et al. 92], 
[Smoot et al. 92]. Today the CMB anisotropy (CMBA) has been detected also 
at intermediate (~ (1° — 10°)) and small angular scales (< 1°), so the CMBA 
angular spectrum is now reasonably known down to the region of the first 
and second Doppler peaks [de Bernardis et al. 00], [de Bernardis et al. 02], 
[Netterfield et al. 02]. Its shape gives information e.g. on the spectrum of 
the primordial cosmological perturbations or can be used to test the in- 
flation theory but rises new questions to which CMBA cannot answers. It 
is however possible to get answers looking at the CMB polarization 
(CMBP). For instance one can use CMBP to disentangle the effects of fun- 
damental cosmological parameters like density of matter, density of dark en- 
ergy etc. This is among the goals of space and ground based experiments 
like [MAP 03], [Delabrouille and Kaplan 02], [Villa et al. 02], [Kesteven 02], 
[Piccirillo et al. 02], [Gervasi et al. 02], [Masi et al. 02] and is the main goal of 
SPOrt a polarization dedicated ASI/ESA space mission on the International 
Space Station [Cortiglioni et al. 99]. 

The relevance of the CMB polarization was remarked for the first time by 
M. Rees [Reees 68]. Since his paper many models of the expected features of 
the CMBP have been published (see for instance [Sazhin 95], [Ng and Ng 96], 
[Melchiorri and Vittorio 97]). They immediately stimulated the search 
for CMBP, but the first detection has been obtained only a few months ago 
([Kovac et al. 02] and [Kogut et al. 03]). Because of its importance this dis- 
covery must be confirmed by new observations made with different systems 
and using different methods of extraction of the CMBP from the sky signal. 
The detection of the CMBP is in fact extremely difficult because the signal is 
at least an order of magnitude below the CMBA level. Moreover polarized fore- 
grounds of galactic origin may cover the CMBP and/or mimic CMBP spots 
by their inhomogeneities. The signal to noise (CMBP / polarized foreground) 
ratio is therefore << 1. In this paper we discuss a method which improves 
this ratio and allows to disentangle CMBP and polarized foregrounds. 

In the microwave range the galactic foregrounds include: 

• synchrotron radiation (strong polarization), 

• free- free emission (null or negligible polarization), 

• dust radiation (polarization possible) 

Beccause here we are interested in polarization, in the following we will ne- 
glect free-free emission. The effects of dust, if present, (e.g. [Prunet et al. 98], 
[Fosalba et al. 02]), will be added to the synchrotron effects. In fact, as it will 
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appears in the following, what matters in our analysis are the statistical prop- 
erties of the spatial distribution of the foregrounds and, by good fortune, the 
spatial distribution of the dust polarized emission is similar to that of the 
synchrotron emission. Behind both types of radiation there is in fact the same 
driving force, the galactic magnetic field which alignes dust grains and guides 
the radiating electrons ([Bernardi et al. 03] and references therein). 

For the anisotropy the separation of foregrounds and CMB was success- 
fully solved by the DMR/COBE team when they discovered CMBA 
[Smoot et al. 92]. The separation of CMBP and foregrounds is more demand- 
ing and recognizing true CMB spots among foreground inhomogneities se- 
vere. Approaches commonly used ( e.g. [Dodelson 97], [Tegmark 99], 
[Stolyarov et al. 02], [Tegmark et al. 96], [Kogut and Hinshaw 00], 
[Bennet et al. 03]), are based on the differences between the frequency spectra 
of foregrounds and CMB, therefore require multifrequency observations. 

In this paper we suggest a different method. It takes advantage of the fact 
that, as we will show in the following, at small angular scales the values of 
the parameters used to describes the polarization of the diffuse radiation at 
a given frequency in different directions fluctuate. We propose of analyzing 
the angular distribution of the polarized radiation on single frequency maps 
of the diffuse radiation and disentangling the main components, polarized 
synchrotron and CMBP, looking to their different statistical properties. This 
method was proposed and briefly discussed in [Sazhin 02]. Here we present a 
more complete analysis. 



2 Polarization Parameters 



2.1 The Stokes Parameters 



Convenient quantities commonly used to describe the polarization 
status of radiation are the Stokes parameters (see for instance 
[Ginzburg and Syrovatskii 65], [Gardner and Whiteoak 66], 
[Ginzburg and Sirovatskii 69]). 

Let's assume a monochromatic plane wave of intensity I and amplitude oc 
\fl. In the observer plane, orthogonal to the direction of propagation of the 
electromagnetic wave, we can choose a pair of orthogonal axes I and r. On 
that plane the amplitude vector of an unpolarized wave moves in a random 
way. On the contrary it describes a figure, the polarization ellipse, when the 
wave is polarized. 
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Fig. 1. The polarization ellipse 



Projecting the wave amplitude on / and f we get two orthogonal, linearly 
polarized, waves of intensity // and I r ,(I = l\ + I r ) whose amplitudes are 
oc \fT\ and oc \[T r respectively. If the original wave of intensity / is polarized, 
Ii and I r are correlated: let's call 1 12 and 7 2 i their correlation products. 



By definition the Stokes parameters are: I = I t + I r , Q = I t — I r , U = Iu + hi, 
and V = i(I 2 i - hi). 

Q and U describe the linear polarization, V the circular polarization and / 
the total intensity. 

The ratio 



tan2% 



U 

Q 



gives the angle x between the vector I and the main axis of the polarization 
ellipse (0 < x < 7r )- 

Rotation by an angle <fi of the (l,r) coordinate system give a new coordinate 
system (l,f) in which the Stokes parameters become 



Q = Q cos 20 — U sin 20 



U = U cos 20 + Q sin 20 



(2) 



So when = x — (1/2) arctan^ the axes of the polarization ellipse coincide 
with the reference axes I and r (see fig. 1). 



4 



2.2 Electric and Magnetic Modes 



To analyze the properties of the CMB polarization it is sometimes convenient 
to use rotationally invariant quantities, like the radiation intensity / and two 
combinations of U and Q: Q + iU and Q — ill. The intensity / can be decom- 
posed into usual (scalar) spherical harmonics Yi m (9,<p). 

I = J2a lm Yi m (9,(p) (3) 

l,m 



The quantities Q ±iU can be decomposed into ±2 spin harmonics 
[Sazhin et al. 96a], [Sazhin et al. 96b], [Seljak and Zaldarriaga 97] 

Q±iu = Y,4£YL 2 {o^) (4) 



The ±2 spin harmonics form a complete orthonormal system (see, for instance, 
[Goldberg et al. 67],[Gelfand et al. 58], [Zerilli 70], [Thorn 80]) and can be writ- 
ten [Sazhin et al. 96b], [Sazhin and Sironi 99]: 

Y l U6,v)=N? m P? m {e)e^ 

(5) 

Y lm 2 (6, V )=N^P lm 2 (e)e^ 



where 

(m + s) (s — m) , . s 

PUx) = (1 " + x) y —^PlT' S - m \x) (6) 



is a generalized Jacobi polynomial, s = ±2 and: 



NL = -1 



'2 1 + 1 (l-s)\(l + s)\ 
4tt \ (l-m)\(l + m)\ 



is a normalization factor. 



The harmonics amplitudes a lm correspond to the Fourier spectrum of the 
angular decomposition of rotationally invariant combinations of Stokes pa- 
rameters. 

4 Alternatively, one can use the equivalent polynomials derived in 
[Goldberg et al. 67] 
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Because spin ±2 spherical functions form a complete orthonormal system: 
J Y± 2 (6, V )Y* n f{6, v )dQ = Uw (7) 



4-7T 



we can write 



af 2 = JdQ (Q(9, <p) ± iU(6, tp)) Y*± 2 (8) 

47T 



Following [Seljak and Zaldarriaga 97] and the very nice introduction made 
more recently by [Zaldarriaga 01] we now introduce the so called E (electric) 
and B (magnetic) modes of these harmonic quantities: 

j£ \ ( -1-2 2\ 

a lm = 2 \ a lm + a lm ) 

(9) 

a fm — 2 — a Zm) 

They have different parities. In fact when we transform the coordinate system 
Oxyz into a new coordinate system Oxyz, such that 



/ = I 



(10) 



the E and B modes transform in a similar way: 



a E = a E 



d B = -a B 



(11) 



Q remains identical in both reference systems and U changes sign. 
It is important to remark that a E and a B are uncorrelated. 
In terms of Q and U we can write: 

°L = \SdSl (Q (Y+ 2 + Y^)+ iU (Y+? - YC m 2 ) ) 
aim = 5/ dSl (iQ (y+ 2 - Yfj) - U (Y+J + Y lm 2 )) 



(12) 
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therefore : 

(KJ 2 ) - ((afj 2 ) 
= 2JdQ((\Q 2 \) - (\U 2 \)) (Y£Y?-* + Y£*Y t -*) 

Here (\Q 2 \) and {\U 2 \) designate values of delta correlated 2D stochastic fields 
Q and U. Omitting mathematical details, the correlation equations for Q and 
U are: 

(QQ*) = \Q 2 \5(n-n') 

(uu*) = \u 2 \5(n - n 1 ) (14) 

(QU*) = 



where d~(Q — Q') = <5(cos# — cos 9') ■ 5(ip — ip') is the Dirac delta - function on 
the sphere. (In the following we will sometimes omit indexes I and m). 



3 Synchrotron Radiation and its Polarization 



Synchrotron radiation results from the helical motion of extremely relativistic 
electrons around the field lines of the galactic magnetic field (see, for instance 
[Ginzburg and Syrovatskii 65], [Ginzburg and Sirovatskii 69], [Westfold 59]). 
The electron angular velocity 

eH p m e c 2 m e c 2 . . 

= 7T- = u —r- (15) 



is determined by the ratio between H p , the component of the magnetic field 
orthogonal to the particle velocity, and £ , the electron energy. As it moves 
around the magnetic field lines the electron radiates. 



b) Single electron 

Until the electron velocity is small (v « c) we speak of cyclotron radiation: 
the electron behaves as a rigid dipole which rotates with gyrofrequency (15) 
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Fig. 2. Spatial distribution of the cyclotron radiation produced by a single electron 




Fig. 3. Velocity cone of an ultrarelativistic (v ~ c) electron (if = magnetic vector, k 
wave vector, tp angle between the electron velocity and the direction of observation) 

in a plane orthogonal to the magnetic field direction and emits a single line. 
The radiation has a dumbell spatial distribution (see fig. 2): 

/(©,$) ~ (1 + cos 2 6), (16) 



is circularly polarized along the dumbell axis (0 = 0) and linearly polarized 
in directions orthogonal to it (0 = 90°). 

When the electron velocity increases the radiation field changes until at v ~ c, 
(£ ^> m e c 2 ) it assumes the peculiar charachters of synchrotron radiation: 

i) radiated power proportional to £ 2 and H 2 , 

ii) continous frequency spectrum, peaked around: 



u m (4>) = u - _ £cQ C s ^ oc H P S 2 (17) 

(ijj is the angle between the velocity vector v and the wave vector k, see fig. 3). 
The peak is so narrow that in a given direction the emission is practically 
monochromatic, 
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Fig. 4. Distribution of the radiation emitted by an ultrarelativistic (v ~ c) electron 
with istantaneous velocity v spiralling (dashed line) around the lines of the magnetic 
field H (orthogonal to the sheet) 

iii)radiation almost entirely emitted in a narrow cone 5 of aperture (see fig. 4) 
$ « — — (18) 



around the forward direction of the electron motion. Inside the cone (cos ip ~ 
1) the frequency is maximum and equal to 

uvo^o(^) , (19) 



In the opposite direction (cos^ ~ —1) intensity and frequency are sharply 
reduced. 

iv)radiation 100 % linearly polarized at the surface of the cone. Inside the cone 
the linear polarization is still dominant but a small fraction of circular polar- 
ization exists, (V ~ o( m ^) and I ~ U ~ Q). Outside the cone the very small 
fraction of radiation produced is elliptically polarized and becomes circularly 
polarized when seen along the H p direction. So the Stokes parameters depend 
on H, the angle fi between H and the line of sight, and the dimensionless 
frequency vjv c , where v c = 1.5(u; m (0)/27r) is the so called critical frequency, 
[Ginzburg and Syrovatskii 65],[Ginzburg and Sirovatskii 69], itself function of 
H (see eq.15). 



b) Cloud of monoenergetic electrons 

When the effects of many monoenergetic electrons with uniform distribution 
of pitch angles are combined, /, Q and U are reinforced (the Stokes parameters 

5 the symmetry plane(@ = 90°) of the cyclotron dumbbell beam, seen by a fast 
moving observer (v f» c) becomes a cone folded around the direction of movement 
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are additive) while V is erased. In fact 
I(u) = Cl H p F(^) 

Q{v) = c 2 H p F p {^) cos2 X 

(20) 

U(v) =c 2 H p F p (%) Bm2 X 
V(u) « 

where C\ and c 2 are constants, F and F p almost monochromatic functions of 
v , and x the angle between the projection of the magnetic vector on the 
observer plane and an axis on that plane (The projection of the magnetic 
vector on the observer plane is the minor axis of the polarization ellipse). If 
inside the emitting cloud the magnetic field varies, also x varies, therefore 
cos 2x and sin 2% must be averaged along the line of sight across the cloud. 
In conclusion the degree of polarization p varies between a maximum value, 
(uniform magnetic field) and zero (magnetic field randomly distributed). 



c) Electrons with power law energy spectrum 

In the interstellar medium the radiating particles are the cosmic ray electrons 
whose energy spectrum : 

N(£) = KE" 1 (21) 

is a power law (see for instance [Gavazzi and Sironi 75], [Longair 94]) and ref- 
erences therein) with spectral index 7 w 2.4 — 3.0 and space density propor- 
tional to K. Because the emission of a single electron is practically monochro- 
matic the resulting radiation spectrum is a power law: 

I(u) = I^H^v-f (22) 

with intensity spectral index (3' = :L ^-, temperature spectral index (3 = (3' + 2 
and Iq a slow function of 7 [Ginzburg and Syrovatskii 65]. (If the magnetic 
field is not uniform we use < H p > instead of H p and a slightly different 
function /o(t))- 
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The Stokes parameters are products of the intensity I{v), the degree of polar- 
ization p and cos 2% or sin 2%, therefore we can write: 



0L±± y-l 

Q = c 3 Hp 2 v 2 cos 2% 



7+1 



7-1 



(23) 



U = c 3 H P 2 



v 2 sin 2% 



V « 



where c 3 is a constant whose value depends on 7 and the distribution of the 
magnetic field along the integration path. 

When the medium inside the synchrotron source, or in the medium where radi- 
ation propagates, is permeated by thermal electrons Faraday rotation modifies 
the polarization charachteristics of the radiation. In fact by Faraday effect, 
when the radiation crosses a region of thickness L permeated by magnetic 
fields and thermal electrons with density n e , the angle of polarization of the 
radiation rotates by an angle (see for instance [Rybicki and Lightman 79]) 



Inside the source this brings depolarization, because radiation produced at 
different points along L suffers different rotations. Additional depolarization 
inside the source comes about when the magnetic field along L is not uniform 
because in eqs.(23)we have to use < sin2x > and < cos 2% > instead of sin2x 
and cos 2%. So the degree of polarization at the source p varies between 
(random magnetic field distribution) and: 



(uniform magnetic field and Faraday rotation absent). 

Outside the source the angle of polarization rotates by Farday effect so the 
Stokes Parameter of the synchrotron radiation Q s and U s measaured by an 
Earth observer are different from the Stokes Parameters at the source. 

A regular trend of the magnetic field inside and outside the source is insuffi- 
cient to guarantee a regular trend of the spatial distribution of Q s and U s . In 
fact as the line of sight moves among adjacent points on the sky n e , < n e > and 



9 = cost < n e > L A 2 



(24) 



Pmax 



3 7 + 3 
3 7 + 7 



< 1 



(25) 
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L fluctuate. So also the angle of Faraday rotation 9 fluctuates by a quantity 



S9 = cost A 2 5(<n e > L) 



(26) 



~ cost A 2 {[L 5(< n e >)] 2 + [<N e > 5(L)] 2 } 1/2 ~ cost A 2 L S(n e ) 



proportional to L. Because L is at least tens or hundred of pc 59 can be 
very large. Besides this effect produced by Faraday rotation, when we look in 
different directions through the interstellar medium Q s and U s may change also 
because H, H p = Hsm/i, 7, f3 and K vary (see for instance [Wielebinski 02], 
[Reich et al. 02]) 6 . 

We conclude that angle and degree of polarization of the synchrotron radia- 
tion fluctuates when the line of sight moves among adjacent regions on the 
sky, a conclusion supported by recent observations which show that the spa- 
tial distribution of the polarized component of the synchrotron radiation at 
intermediate and small angular scales is highly structured. On maps we see in 
fact lines (channels) along which p goes to zero and sudden rotations of the 
plane of polarization when the line of sight goes from a side to the other of a 
channel [Shukurov and Berkhuijsen 03]. 

We expect therefore that at small and intermediate angular scales Q s and U s 
are stochastic functions of the direction of observations so we can write: 



Exceptions to this behaviour can be expected when along the line of sight 
there are peculiar regions characterized by special field configurations and/or 
by large densities of thermal electrons, like in the well known loops and spurs 
or in HII regions, where systematic effects overcome random effects. Nature 
and origin of these features, clearly visible on Brouw and Spoelstra maps of 
the polarized component of the galactic background [Brouw and Spoelstra 76], 
are discussed for instance by [Salter and Brown 88]. Usually these regions are 
close to the observer and have large angular dimensions, so their effects can be 
removed if one evaluates the angular power spectrum of the radiation distri- 
bution and limits his analysis to angular scales below few degrees (multipole 
order / > 150) 7 . 

6 e.g. F = (H sin u)^~ is a nonzero mean, nonzero variance variable; cos 2% and 
sin 2x are zero mean, nonzero variance variables 

7 At larger angular scales (smaller values of 1) the number of independent samples 
one collects decreases until it is insufficient to carry on statistical analysis and the 
method of separation of CMB and foregrounds, we will illustrate in the following, 




(27) 
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Being aware that this is a very important issue we tested its validity analyzing 
the distribution of the measured values of Q and U on maps extracted from 
the Effelsberg survey ([MIPFR-Bonn 03] and [Duncan et al. 99]) of the diffuse 
radio emission. We found (see Appendix 2) that up to angular scales of 5 
degrees (1 close to 36) the difference < Q 2 > — < U 2 > is fully consistent with 
zero. 

Because of eq. (27) if Q s ^ also U s ^ 0, so the synchrotron radiation has 
both electric and magnetic modes: 

a E ' s ± 0, a B ' s ^ (28) 

and (see eqs. (13) and (27)): 

((a E n 2 ) = ((a B n 2 ) (29) 



4 Stokes Parameters of the CMB 



In a homogeneous and isotropic Universe only temperature and intensity 
I — Ii + I r change as the Universe expands : both decrease adiabatically. 
Because this is true for I t and I r separately, we do not expect anisotropy nor 
polarization therefore Q = and U — I u — are natural consequences. 

On the contrary, inhomogeneities and perturbations of matter density or of 
gravitational field, induce anisotropy and polarization of the CMB. At the re- 
combination epoch linear polarization appears as a by product of the Thom- 
son scattering of the CMB on the free electrons of the primordial plasma. The 
polarizarion tensor it gives can be calculated solving the Boltzman transfer 
equations of the radiation in a nonstationary plasma permeated by a variable 
and inhomogeneous gravitational field [Basko and Polnarev 80], [Sazhin 84], 
[Harrari and Zaldarriaga 93], [Sazhin 95]. 

The gravitational field is made of a background field, with homogeneous and 
isotropic FRW metric, and an inhomogeneous and variable mix of waves: den- 
sity fluctuations, velocity fluctuations, and gravitational waves. Because of 
their transformation laws these waves are also said scalar, vector and tensor 
perturbations, respectively. 



fails. 
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4-1 Scalar (density) perturbations 

Scalar (density) perturbations affect the gravitational field, the density of 
matter and its velocity distribution. They were discovered studying the matter 
distribution in our Universe on scales from ~ 1 Mpc to ~ 100 Mpc. It is firmly 
believed they are the seeds of the large scale structure of the Universe and are 
reflected by the large scale CMB anisotropy detected for the first time at the 
beginning of the '90s [Strukov et al. 92], [Smoot et al. 92]. Their existence is 
predicted by the great majority of models of the early Universe. 

Observation shows that the effects of these perturbations are small, so we 
can treat them as small variations 5i, 5 r , 5 U of I r and /„. Introducing the 
auxiliary functions a and (3: 

5i + 5 r = (/i 2 - 

5 l -5 r = (l- f i 2 )(3, (30) 

(/i is the angle between the line of sight and the wave vector) for plane 
waves the Boltzman equations (see for instance [Sazhin 95], [Sazhin et al. 96a], 
[Sazhin et al. 96b] and reference therein) become: 

% = F - ^a T n e a(r])a - ^a T n e a(r])p 

(31) 

% = -^cr T n e a(rj)a - ^a T n e a(r])(3 

where F is the gravitational force which drives both anisotropy and polariza- 
tion, o T is the Thomson cross-section, n e is the density of free electrons, and 
a{r\) the scale factor. 

These equations give: 

Q = -l(l-^)JF(r ] ) (e-r-e-fr)dri 

(32) 

U = 
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where t{j]) is the optical depth of the region where the phenomenon occurs. 
Rotating the coordinate system we can generate a new pair of Stokes param- 
eters (Q',U r ): no matter which is the system of reference we choose these 
parameters satisfy the symmetry parity condition. 

Because there is always a system in which Q ^ and U = 0, we may conclude 
that, in a system dominated by the primordial density perturbations (see, for 
instance [Dolgov et al. 90]) magnetic modes of the CMB polarization vanish 
and only electric modes exist [Seljak and Zaldarriaga 97]. Therefore: 

a E > d ^ a B4 = (33) 



where index d stays for density perturbation. 



4-2 Vector (velocity) perturbations 

Vector perturbations, associated to rotational effects, perturb only velocity 
and gravitational field. They are not predicted by the inflation theory and it is 
common believe that they do not contribute to the anisotropy and polarization 
of the CMB. 



4-3 Tensor (gravitational waves) perturbations 

Gravitational waves (tensor modes) and gravitational lensing of large scale 
structures of the Universe induce B modes in the distribution of CMBP(see 
for instance [Hu et al. 03])) . Gravitational lensing gives a power spectrum at 
least an order of magnitude below the power spectrum of the E mode polarized 
signal produced by scalar perturbations, with a maximum at I ~ 1000). The 
power spectrum of the B mode polarization produced by gravity waves is 
maximun at I ~ 90, is definitely below the power spectrum produced by 
gravitational lensing for / > 100 and only at very large angular scales (I ~ 20) 
it may be comparable to the scalar E modes. 



If one excludes very large angular scales we may conclude that CMBP is 
dominated by E-modes. B-modes are just a contamination by B-modes at 
levels of 10% or less. 
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5 Separation of the polarized components of Synchrotron and CMB 
Radiation 

The CMB radiation we receive is mixed with foregrounds of local origin. When 
the anisotropy of the CMB was detected, to remove the foregrounds from maps 
of the diffuse radiation, data were reorganized in the following way: 

f d (x, y) = fjx, y) + Mi(x, y)f ijC (x, y) (34) 

where T^, is a two dimension vector (map) which gives the total signal mea- 
sured at different points (x, y) on the sky, T n the noise vector, M the matrix 
which combines the components T itC of the signal. At each point (x, y) we can 
in fact write: 

Td QrxXn 9cmbT C mb,c Qsyn^synfi Qff^ffjC 9dustJ^dust,c ••• (^5) 

where gi are weights, given by M, T synfi is the synchrotron component, 2// )C 
the free- free emission component, Td us t,c the dust contribution and so on. Using 
just one map the signal components cannot be disentangled. If however one 
has maps of the same region of sky made at different frequencies it is possible 
to write a system of equations. Provided the number of maps and equations 
is sufficient, the system can be solved and the components of T d separated, 
breaking the degeneracy. We end up with a map of T cm b which can be used to 
estimate the CMB anisotropy. 

When we look for polarization at each point on the sky we measure tensors 
instead of scalar quantities, therefore to disentangle the polarized components 
of the CMB we need a greater number of equations. Here we will concentrate 
on the separation of the two dominant components of the polarized diffuse 
radiation: galactic synchrotron (plus dust) foreground and CMBP. 



5. 1 The estimator D 

Instead of observing the same region of sky at many frequencies, we suggest a 
different approach. It takes advantage of the differences between the statistical 
properties of the two most important components of the polarized diffuse 
radiation (CMB (background) and synchrotron (foreground) radiation) and 
does not require multifrequency maps. 
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We define the estimator: 

^((«T)-((aT) (36) 



where a E = a E,s + a E,c and a B = a ByS + a B ' c (here and in the following indexes 
s or c stay for synchrotron and CMB, respectively). 

Because 

• E and B modes of synchrotron do not correlate each other neither correlate 
with the CMB modes 

((a E ) 2 ) = ((a E > s ) 2 ) + ((a E > c ) 2 ) + 2 (a E ' s a E > c ) = ((a E > s ) 2 ) + ((a E > c ) 2 ), 

((a B ) 2 ) = ((a B > s ) 2 ) + ((a B < c ) 2 ) + 2 (a B > s a B > c ) = ((a B > s ) 2 ) + ((a B ' c ) 2 ), 

• < a E ' s >=< a B > s >^ 

a E, c = a E,d^ Q a B ' c = a B ^ < 0.1 a E > d 

where indexes d and t stay for density (or scalar) and tensor perturbation, 
respectively. 



D gives an estimate of the E-mode excess in maps of the polarized diffuse radi- 
ation. If tensor perturbations are negligible this excess is the CMBP signal. If 
tensor perturbations are important the excess is a lower limit with a system- 
atic difference from the true value which in the worst condition (maximum 
contribution to CMBP from gravitational waves and gravitational lensing) 
reaches a maximum value of 10%. 

Let' s now consider the angular power spectrum of D. For multipole / we can 
write: 

A = (a?) 2 - {aff = -±- £ (\a E J 2 - \a B m \ 2 ) (37) 

' rn=—l 



where \af m \ 2 and \af m \ 2 are random variables with gaussian distribution p(af^f) 
(see eqs.(56) and (57) in Appendix A). According the ergodic theorem (in 
the limit of infinite maps, the average over 2D space is equivalent to the 
average over realisations) the average value of Di is equal to the difference of 
the average values of |o^J 2 and |o^J 2 summed over m. Taking into account 
equation (57) we can therefore write: 

(A) = (\af | 2 - |af| 2 ) (38) 
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where 8 : 

oo , 

(«f) 2 = / WF m \ 2 P(aFJd 2 af m = (\af m \ 2 



(39) 



oo 

B\2 _ r \„B \2^(„B \J2„B _ \„B 12 



K) 2 = / 



lm — \nrri 



Comparing eq.(39) with the ordinary definition of multipole coefficients: 



1 1 

C E,B — ^ /\„ E < B \ 2 



21 + , 

m=—l 



T E (K m T) (40) 

-L , 7 



we can write: 



(af f = Cf 



if? 



(41) 



5.2 Separation uncertainty 



For synchrotron radiation (Di) should be zero, non zero for density pertur- 
bations, but on real maps it is always different from zero. In fact a map is 
just a realization of a stochastic process and the amplitudes of \af m \ 2 and 
|a^| 2 , averaged over 2D sphere, have uncertainties which add quadratically, 
so (Di) 7^ even in the case of synchrotron polarization. This effect, very sim- 
ilar to the well known cosmic variance of anisotropy [Knox and Turner 94], 
[Sazhin et al. 95] (the real Universe is just a realisation of a stochastic process, 
therefore there will be always a difference between the realization we measure 
and the expectation value) does not vanish if observations are repeated. 



The variance of D t is 

v(A 2 ) = (A 2 ) - (A) 2 (42) 

The quantities a E and a B , being sums over m of 21 + 1 stochastic values with 

8 Equation (39) is an explicit form of the average of the stochastic variables |a^,| 2 
and \af m \ 2 over a probability density p(a lr ^ ), the short form being triangle brackets 
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gaussian distribution, have a x 2 distribution with 2*(2/+l) degrees of freedom, 
so their variance is 

s ^ « i/irr (43) 

More explicitly 

(A 2 ) = 2[^Y ((«f ) 4 + («f ) 4 ) ( 44 ) 
When the synchrotron foreground is dominant (af = af ) 

in agreement with (43). 



5. 5 A criterium for CMBP detection 



The synchrotron foreground is a sort of system noise which hampers the de- 
tection of the signal, the CMB polarization. At frequencies sufficiently high 
(above ~ 50 GHz, see next section) the noise is small compared to the sig- 
nal therefore direct detection of CMBP is possible. At low frequencies on the 
contrary the CMBP signal is buried in the noise created by synchrotron and 
dust emission. In this case to recognize the presence of the CMB polarization 
we can use our estimator D. 

At angular scale /, to be detectable the CMBP must satisfy the condition 

Cf ' c > A ■ Cf' s (46) 



where Cf' 1 are the coefficients of the multipole expansion of the E modes and 
A is the confindence level of the signal detectability. In a similar way we can 
write for our estimator: 

Df ' c > A ■ Df' s (47) 



where 



^ = V2TTI (af ' S)2 (48) 
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If one neglects tensor perturbations (Df ' c = Cf' d ) the criterium for the CMBP 
detection becomes 




D?*>A.J-—C*>> (49) 



so by D we get the CMBP level with an uncertainty ~ \J Wp[C^ s which 
decreases as I and the angular resolution increase. Tensor perturbations, if 
present, add to this uncertainty a systematic uncertainty < 10%. 



6 Angular power spectra of polarized synchrotron. 



The angular power spectra of the polarized component of the synchrotron radi- 
ation have been studied by [Baccigalupi et al. 02], [Burigana and La Porta 02], 
[Tucci et al. 02a], [Tucci et al. 02b], [Bruscoli et al. 02] and [Giardino et al. 02] 
using the few partial maps of the polarized diffuse radiation available in lit- 
erature. It appears that the power spectra of the degree of polarization p and 
of electric and magnetic modes E and B follow power laws of I up to I ~ 10 3 . 
For the degree of polarization the spectral index is a p ~ 1.6 — 1.8. For the E 
and B modes different authors get different values of the spectral index. 

According the authors of papers [Tucci et al. 02b] and [Bruscoli et al. 02], 
Parkes data give: 

i \ a E .B 

(50) 




with cxe ~ olb ~ 1.4 -r- 1.5 and dependence of a on the region of sky and the 
frequency. 

In paper [Baccigalupi et al. 02], using Effelsberg and Parkes data, the authors 
get: 

Cf^C .10-(^)".(^) 2S (51) 



with a = (1.8±0.3), (3 = 2.9 and C = (1.6 ± 1) (here we adjusted the original 
expression given in [Baccigalupi et al. 02] writing it in adimensional form). 

In paper [Giardino et al. 02] the authors, using a completely different set of 
observational data ([Jonas et al. 98] and [Giardino et al. 01]) conclude that in 
the multipole range / = 40—250 the spectral indexes of the E and B modes are 
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dE = (1.57 ±012) and cub = (1.45 ±0.12) respectively, while for the polarized 
intensity the spectral index is aj = (2.37 ± 0.21). 

Inside the multipole range (10 2 < I < 10 3 ) the values of ole,b obtained by 
the three groups are marginally consistent and reasonably close to ~ 1.8. 
At lower value of I definite differences exist between the results obtained by 
[Baccigalupi et al. 02] and by [Giardino et al. 02]: these differences must be 
understood, but they are not important here because our statistical method 
cannot be used for small values of / when the number of samples becomes in- 
sufficient to carry on statistical analyses. For high values of / , if one excludes 
extrapolations and models (eg. [Giardino et al. 02]), there are no data in lit- 
erature, but this is not a limitation because for large values of / the results of 
our method are practically independent from the spectral index. 

All the above authors got their results analyzing low frequency data (1.4, 2.4 
and 2.7 GHz), therefore the extension of their spectra to tens of GHz, the 
region where CMB observations are usually made, depend on the accuracy 
of P, the temperature spectral index of the galactic synchrotron radiation. 
A common choice is (3 = 2.9 but in literature there are values of (3 rang- 
ing between 2.5 and 3.5. Moreover f3 depends on the frequency and 
the region of sky where measurements are made (see [Gavazzi and Sironi 75], 
[Salter and Brown 88], [Zannoni et al. 00], [Platania et al. 98]). Last but not 
least many of the values of (3 in literature have been obtained measuring the 
total (polarized plus unpolarized) galactic emission. In absence of Faraday 
effect (3,/3 po i and (3 unpo i, the spectral indexes of the total, polarized and un- 
polarized components of the galactic emission, should be identical (see eqs. 
(22), (23)). However when Faraday effect with its v~ 2 frequency dependence 
is present, we expect that, as frequency increases, the measured value of the 
degree of polarization (se eq. (25)) increases, up to 

_ 3(3-3 

V — Pmax „ a i 



therefore we should observe (3 po i < (3. The expected differences are however 
well inside the error bars of the data in literature so at present we can neglect 
them and assume (3 ~ f3 poi ~ (3 unpo i. 

Instead of extrapolating low frequency results it would be better to look for 
direct observations of the galactic emission and its polarized component at 
higher frequencies. Unfortunately above 5 GHz observations of the galactic 
synchrotron spectrum and its distribution are rare and incomplete. At 33 GHz 
observations by [Davies and Wilkinson 99] give at some patches on the sky a 
galactic temperature of about 2 AjiK from which follows that at the same 
frequency we can expect polarized foreground signals up to several \iK. At 
14.5 GHz observations made at OVRO [Mukherjee et al. 03] give synchrotron 
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Fig. 5. Power spectrum versus multipole order I of the polarized components of CMB 
(solid line) and three possible power spectra of the galactic synchrotron radiation 
(q = 1.5 (dotted line), a = 1.8 (dashed line) and a = 2.1 (dot-dash line)) calculated 
at 43 GHz. (See text for details of calculations and model) 

signals of 175 fiK, equivalent to about 15 fiK at 33 GHz of which up to 10 
\iK can be polarized. 

In conclusion there are large uncertainties on the frequency above which ob- 
servations of the CMBP are practically unaffected by the polarized compo- 
nent of the galactic diffuse emission. To be on the safe side we can set it at 
50 GHz (see for instance [Baccigalupi et al. 02], [Burigana and La Porta 02], 
[Tucci et al. 02a], [Tucci et al. 02b] and [Bruscoli et al. 02]). Above 50 GHz 
CMBP definitely overcomes the polarized synchrotron foreground. Below 50 
GHz contamination by the galactic emission can be important but its evalua- 
tion, usually made by multifrequency observations, is dubious. 

This conclusion is supported by figure 5 where we plotted, versus the multipole 
order /, the power spectra of CMBP and galactic synchrotron at 43 GHz. The 
CMB spectrum has been calculated by CMBFAST [Seljak and Zaldarriaga 98] 
assuming standard cosmological conditions (CMB power spectrum normalized 
to the COBE data at low I, £l b = 0.05, Q C dm = 0.3 , fi A = 0.65, Vt v = 0, H = 
65 km/sec/Mpc, Tqmb = 2.726K , Yn e = 0.24, standard recombination). The 
synchrotron spectrum has been calculated assuming for E\ and Bi the scaling 
law (51) with C = 2.6 (most pessimistic case), a = 1.5; 1.8; 2.1 and (3 = 2.9 
respectively. It appears that at 43 GHz the CMB power is comparable to the 
synchrotron power only at very small angular scales (I > 500). 

Similar calculations at other frequencies confirm that only above ~ 50 GHz 
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and at small angular scales (large values of I) the CMBP power spectrum 
overcomes the synchrotron spectrum. Below ~ 50 GHz direct detection of the 
CMBP is almost impossible even at small angular scale. 



To overcome this limit we can use our estimator. To see how it improves the 
CMBP detectability we calculated the angular power spectrum of Df, the 
estimator we expect when the diffuse radiation is dominated by synchrotron 
radiation. From eqs.(40), (45) and (51) we get: 



a quantity which can be directly compared with the power spectrum of CMBP. 
In fact eqs.(33), (38), (41) and (44) show that the power spectrum Df of the 
estimator evaluated when the sky is dominated by CMBP density pertur- 
bations coincides with the power spectrum of CMBP produced by density 
perturbations. 

Figure 6 shows at 37 GHz the power spectra at 37 GHz of: i)Df, the estimator 
for a sky dominated by galactic synchrotron (solid line , calculated using 
eq.(51)), n)Df, the estimator for a sky dominated by CMBP. It coincides 
with the power spectrum of CMBP (dotted line, calculated as in figure 5 with 
CMBFAST using the same standard cosmological conditions), iii)the power 
spectrum of the polarized component of the galactic synchrotron radiation 
(dashed line, calculated using eq.(51)). As expected at 37 GHz the CMBP 
power is well below the synchrotron power, therefore direct observations of 
CMBP are impossible (the maximum value of the CMBP power is about 2.5 
times below the synchrotron power at the same /). However above / ~ 250 the 
power of the synchrotron estimator is definitely below the power of the CMBP 
estimator: at I = 1000 the ratio CMBP/D is maximum and close to 7. This 
confirms that the use of D allows to recognize the CMBP also at frequencies 
well below 50 GHz. 



7 Simulations 



To further test the capability of our estimator we studied the separation of 
CMB and galactic synchrotron using measured instead of expected values 
of Z^.The measurements were simulated by random numbers, with gaussian 
distribution ~ JV(0, 1), zero mean and unity variance (see eqs (56) - (57)). 
Two series of 21 + 1 random numbers gave representations of af m and af m 
respectively and from them we obtained D\ (see eq.(37)). 
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Fig. 6. Power spectrum at 37 GHz of the expected value of the estimator for a 
galactic synchrotron dominated sky (solid line) . The dotted line represents both the 
estimator and the CMBP power spectrum for a CMBP dominated sky. The dashed 
line is the power spectrum of the polarized component of the galactic synchrotron 
radiation (see text for details of model and calculations) 
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Fig. 7. Multipole power spectrum of simulated measurements of the estimator D 
for a synchrotron dominated sky (see text) with infinite angular resolution (no 
smoothing on I, (A/ = 1)) 
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Fig. 8. Same as figure 7 with finite angular resolution (smoothing on Al = 10) 

Then, to take into account that real data are collected with a finite angu- 
lar resolution (e.g. Boomerang data come from regions whose angular exten- 
sion is equivalent to Al ~ 100 [de Bernardis et al. 00], [de Bernardis et al. 02], 
[Netterfield et al. 02]) we averaged the above values of A on intervals Al = n: 

k=+n 

Di = — rl E A-n/2+Jkl (53) 

U 1 k=0 



Figure 7, figure 8 and figure 9 show \Di\ (the sign of D\ is arbitrary) versus I, 
for Al = 1, Al = 10 and Al = 100 respectively: the very large fluctuations of 
the estimator are drastically reduced as soon as Al increases. 

Figure 10, figure 11 and figure 12 are similar to figure 6. Here, instead of the 
expectation value, we plot simulated measurements of Df at 37 GHz, for Al = 
1 (figure 10), Al = 10 (figure 11) and Al = 100 (figure 12), respectively. Once 
again the CMBP power spectrum comes from CMBFAST assuming the same 
cosmological conditions we assumed above. The synchrotron power spectrum 
is obtained from eq.(51) with C Q = 2.6, a = 1.5,(3 = 2.9 (most pessimistic 
condition). 

Figure 13 shows at 17 GHz the same quantities we plotted in figure 12. For a 
better appreciation of the differences among the three curves, on the vertical 
axis here we use a logarithmic scale. The power spectrum of the estimator 
now almost touches the two highest peaks of the CMBP spectrum. Probably 
17 GHZ is the lowest frequency at which, in the most favorable conditions, 
one can use D t . In the most pessimistic case (C = 2.6 and a = 1.5) the 
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Fig. 9. Same as figure 7 with finite angular resolution (smoothing on Al = 100) 
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Fig. 10. Similar to figure 6. Here we plot simulated measurements with infinite angu- 
lar resolution (no smoothing on I) instead of the expectation value of the estimator 
at 37 GHz for a synchrotron dominated sky. The dotted line give both the estima- 
tor and the CMBP power spectrum for a CMBP dominated sky. The dashed line 
is the power spectrum of the galacic synchrotron when its expected contribution is 
maximum (eq. (51) with j3 = 2.9, a = 1.5 and C Q = 2.6) (see text for details) 



corresponding frequency is 25 GHz. 
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Fig. 11. Same as figure 10 with finite angular resolution (smoothing on Al = 10) 
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Fig. 12. Same as figure 10 with finite angular resolution (smoothing on Al = 100) 
8 Impact of real world experimental conditions 



We analyzed the impact of the real world experimental conditions on our 
method for disentangling CMBP and foregrounds in maps of the polarized 
diffuse radiation. 

A polarimeter is a two channel system which (by hardware and/or software 
methods) splits the sky signal in two polarized components, send them to 
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Fig. 13. Same as figure 12 at 17 GHz using logarithmic scale on the vertical axis 

separate channels where they are amplified and then looks for correlations 
between the two components. The system outputs are proportional to a pair 
of Stokes parameters (e.g. U and Q), or to the electric and magnetic modes a E 
and a B or to a combination of them, e.g. our estimator D. Unfortunately as the 
signals propagate through the system, noise is added to them so at the system 
output the signal is mixed to noise. Moreover if the channels are asymmetric 
or there are cross talks between them, the outputs contain additional signals 
which simulate spurious polarization, usually an offset of the system output 
from the level one should expect with an ideal system when polarization is 
absent (for a discussion see for instance [Spiga et al. 02]). So in the real world 
the signal we are looking for, S trU e, is accompanied by uncorrected noise and 
offset. What one measures is 

Smeas = Strue + V /=== ~ + ^(S trU e + N) (54) 



where 77 ~ 1, e is the cross talk coefficient, N the noise signal and a noise the 
noise standard deviation. In commercial systems e ~ (1CT 2 — 1CT 3 ) while 
in dedicated CMBP experiments values of e < 10 -5 have been obtained 
[Cortiglioni et al. 99]. By accurate choice of the system components we can 
therefore minimize the offset and further reduce it by phase modulation tech- 
niques (see e.g. [Spiga et al. 02]). 

At this point we have noisy maps we can use to evaluate D and its power 
spectrum instead of U and Q or a E and a B . As shown by the above simulation 
the signal/noise ratio for D is definitely better than the signal/noise ratio 
for the Stokes Parameters or a E or a b and brings it to values similar to the 
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ones we find when studying weak radiosources or, in the worst situation, the 
CMB anisotropy. From here we can therefore go on and extract D from the 
remaining noise using the well known methods of time integration, commonly 
used in radioastronomy. 



9 Conclusions 

Observations of the CMB polarization are hampered by the presence of a 
foreground, the polarized component of the galactic synchrotron radiation. 
Only above ~ 50GHz the cosmic signal definitely overcomes the galactic syn- 
chrotron and direct measurements of the CMBP are possible. Between ~ 30 
GHz and ~ 50 GHz background and foreground are comparable. Below ~ 30 
GHz the polarized sky is dominated by the galactic signal. 

So when measurements are made at ground observatories where atmospheric 
absorption prevent observations above ~ 40 GHz, all programs for CMBP 
measurements must include observing and analysis strategies for disentangling 
the galactic synchrotron signals from the CMBP signals. A common approach 
is fitting models of the intensity, frequency dependence and spatial distribution 
of the cosmic and galactic signals to multifrequency maps of the polarized 
diffuse radiation. Or, when observation are made at one frequency, looking for 
additional data in literature, but here the accuracy of the available data is 
insufficient to get firm evaluations of CMBP. 

In this paper we presented a different approach which at small angular scales 
(< 0.7° (/ > 250)) and down to frequencies as low as ~ 25GHz (~ YlGHz in 
the most favorable conditions) allows to extract the CMBP signal from single 
frequency maps of the polarized diffuse radiation. It takes advantage of the 
different statistical properties of the spatial distributions of CMBP and polar- 
ized galactic synchrotron. By our estimator D, which evaluates the difference 
between E- and B-modes, we get the polarized component of CMB with a 
maximum systematic (underestimate) uncertainty of 10%. This uncertainty is 
set by the contamination by the tensor perturbations which add B-modes to a 
CMBP map dominated by the E-modes generated by scalar(density) pertur- 
bations. Improving our knowledge of the tensor perturbations we will reduce 
the above uncertainty and improve the accuracy of our method of measuring 
CMBP. 

The accuracy we can get with our method is the maximun one can obtain 
at ground observatories with today (2nd generation) systems for measuring 
the CMB polarization. These 2nd generation experiments are just arrived on 
the verge of detecting the polarized signals produced by density perturbations 
(see for instance [Kovac et al. 02] and [Kogut et al. 03]). Direct observations 
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of the signals associated to tensor perturbations requires new, 3rd generation, 
intrinsically able to reject the foreground signals which, at present, are still in 
preparation (see for instance [Piccirillo and Giraud-Heraud 03]). 
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Appendix 1: Stochastic properties of the harmonics amplitudes 

Here and overall in paper we suppose that af m are complex random variables 
which satisfy the probability distribution law: 

K«fJ = d? exp (-*§£) 

(55) 

p(afj = ^exp 

I \ Im / 

with variance (\af m \ 2 ) = E? and (\af m \ 2 ) = Bf. 

They have all the propertiers of gaussian variables (below we omitt indexes E 
and B in first and second equations): 

oo 

/ p(aim)d 2 ai m = 1 

— oo 

(56) 

oo 

/ aimP{aim)d 2 ai m = 

— oo 

7 \af m \MaFJd 2 af m = E 2 

— oo 

(57) 

7 \aL\MaL)d 2 at = Bf 

— oo 
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Setting 

oo 

(F) = J Fp(a lm )d 2 a lm (58) 

— oo 

with current index E or B, it immediately follows: 

■ oo 

lafj 4 ) = / \a?MafJd 2 af m = 2Ef 

— oo 

(59) 

KJ 4 ) = 1 \«LMo»J*af n = 2Bf 



Appendix 2: Stochastic properties of the synchrotron radiation at 
low galactic latitudes 



In the our paper we assert that for the galactic synchrotron emission the 
measured values of the Stokes Parameters Q and U behave as stochastic vari- 
ables and random fields with gaussian distribution. This statement certainly 
holds at high galactic latitudes. At low galactic latitudes, where we observe 
large scale galactic structures, regular magnetic fields [Duncan et al. 99] and 
quasi-periodic structures with typical sizes of about 250 pc, (the amplitudes 
of regular and irregular components are approximately equal) this assumption 
has to be checked. 

To do it we analyzed the distribution of the measured values of Q and U 
on regions of different extension extracted from the Effelsberg maps of the 
polarized diffuse radiation [MIPFR-Bonn 03]. The six fields were chosen within 
±5° from the galactic plane, at galactic longitudes between 68° and 16°. They 
had dimensions of 0.5° x 0.5°, 1° x 1°, 2° x 2°, 3° x 3°, 4° x 4°, 5° x 5°, 
respectively 

For each region we examined the shapes of the distributions, and calculated 
average value and the variance of the measured values of Q and U. The dis- 
tribution of Q in a 5° x 5° area is shown in Figure 14. Very similar results are 
obtained for U and for smaller regions. When the region is very small, (Figure 
15 is for an 0.5° x 0.5° area), the distribution becomes broad, but its shape is 
still gaussian. Results of analyses carried on on all the areas are presented in 
Table 1: for sufficiently large samples both Q and U are compatible with zero. 
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Table 1 

Characteristics of the observed distributions of Q and U measured at low galactic 
latitudes in areas of various extensions (the larger is the number of pixels the closer 
is the histogramm to a gaussian shape) 



size 


average 


variance 


average 


variance 


total 


of 


value 


of 


value 


of 


number 


patch 


of Q 


Q 


of U 


U 


of pixels 


0.5° x 0.5° 


16.2 


564.0 


16.9 


733.4 


256 


1° x 1° 


6.2 


269.3 


-6.2 


311.9 


961 


2° x 2° 


0.07 


426.7 


0.17 


478.7 


3721 


3° x 3° 


2.3 


675.0 


-3.2 


498.7 


8281 


4° x 4° 


-0.95 


372.0 


-2.1 


357.6 


14641 


5° x 5° 


-1.4 


528.7 


-1.8 


507.7 
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Fig. 14. Distribution of the values of the Stokes Parameter Q of the diffuse radiation 
measured in a 5° x 5° area of sky within ±5° from the galactic plane 
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Fig. 15. Distribution of the values of the Stokes Parameter Q of the diffuse radiation 
measured in a 0.5° x 0.5° area of sky within ±5° from the galactic plane 



To check the validity of eqs. (42) to (45) we evaluated the variance 



5D 2 = j2/Na 2 QtU 



(60) 



where N is the number of pixels in the field, Oq v is the measured variance of 
Q and U and D ~ (Q 2 ) — (U 2 ) is our estimator. The quantity (60) is plotted 
in Figure 16 versus the number N of pixels. In the same figure we plot also: 
i)the calculated values of the differences (Q 2 ) — (U 2 ) in each field and ii)the 
calculated values of 



fi/N 



((Q - (Q)) 2 ) - ((u - (U)) 2 } 



(61) 
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Fig. 16. Plot versus the number N of pixels of the expected (dashed line, eq.(60), a 
estimated as sum of the variances of Q and U on all the studied fields) and measured 
(solid line, eq. (61)) values of the variance of our estimator T)D. Diamonds are the 
calculated differences between variances of Q and U and actual data. (Logarithmic 
scale are used for clarity). 
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